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Using numerical techniques we study the spectral function A(k, lu) of a spin-fermion model for 
cuprates in the regime where magnetic and charge domains (stripes) are developed upon hole- 
doping. /.From A(k,Lj) we study the electronic dynamics and determine the Fermi Surface (FS), 
which is compared with angular resolved photoemission results for La2-xSrxCu02. A pseudogap 
is observed in the density of states at the chemical potential for all finite dopings. The striped 
ground state appears to be metallic in this model since there is finite spectral weight at the chemical 
potential, but the electronic hopping seems to be stronger perpendicular to the stripes rather than 
along them. The band structure is not rigid, contrary to the behavior found in mean-field studies, 
and changes with doping. Both mid-gap (stripe induced) and valence band states determine the FS. 
For vertical (horizontal) stripes, a clear FS appears close to (7r,0) ((0, vr)), while no FS is observed 
close to (0, tt) ((tt, 0)). Along the diagonal direction the spectral function shows a clear quasi-particle 
peak close to (0,0), but its weight is reduced as the chemical potential is approached. A weak FS 
develops along this direction as the system is doped. 

PACS numbers; 71.10.Fd, 74.20.Mn, 74.72.-h 



I. INTRODUCTION 

Neutron scattering studies have shown that many high 
temperature superconducting materials exhibit magnetic 
incommensurabihty (IC) upon hole doping. |l] In some 
compounds it is believed that this spin IC is due to 
the formation of charge stripes, j^] Indirect evidence of 
static stripes formation has been observed in Nd doped 
Laz-xSr^CuOs (LSCO). § On the other hand, in LSCO, 
angular resolved photoemission experiments (ARPES) 
1^ as well as neutron scattering , show features consis- 
tent with the existence of dynamical stripes. In neutron 
scattering the evidence is in the broadening of the peak 
in the supposedly inhomogeneous stripe regime, while 
in ARPES a one-dimensional-like band dispersion and 
a straight Fermi surface (FS) around (tt, 0) and (0, tt) are 
observed. 

The goal of this paper is to gain theoretical under- 
standing on the effect of stripes in the spectral func- 
tion and FS of a system of electrons and spins, since 
these are properties which can be measured in ARPES 
experiments. In previous attempts to study these prop- 
erties, mean-field and exact diagonalization (ED) tech- 
niques have been applied to the Hubbard and t-J mod- 
els. The problem with these models is that it has 
not been shown that their actual ground-state has striped 
characteristics. In mean field approaches § not all possi- 
ble states have been considered, and with ED, additional 
attractive terms have to be added to the Hamiltonian to 
stabilize the striped phase, p] Thus, here we consider a 



spin-fermion model which can be unbiasedly studied with 
accurate Monte Carlo (MC) techniques without "sign 
problems" and which, upon hole-doping, has a striped 
ground-state, as shown in previous investigations. 
This model has been used successfully to qualitatively 
describe magnetic and charge properties of the cuprates, 
showing that spin IC and charge stripe formation are re- 
lated. In 

The main result of the present effort is that mid-gap 
states, associated with the stripes, and valence band 
states related to the background, both contribute to de- 
termine the FS. The band structure is not rigid, as found 
in mean-field calculations. [D A pseudogap at the chem- 
ical potential clearly develops as the system is doped. 
The doped holes are not only introduced into the stripes 
but also into the background, and the striped state ap- 
pears to be metallic, although the electronic hopping is 
stronger along the direction perpendicular to the stripes. 

The paper is organized as follows: the spin-fermion 
model is presented in section II, an schematic picture is 
shown in section III, and section IV is devoted to the 
results. The pseudogap is discussed in section V, and 
the conclusions are presented in section VI. 

II. THE MODEL 

The spin-fermion model is constructed as an interact- 
ing system of electrons and spins, crudely mimicking phe- 
nomenologically the coexistence of charge and spin de- 
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grees of freedom in the cuprates. [p[-pl|. Its Hamiltonian 
is given by 



h.c.) + J^Si-Si 



(ij> 



(1) 

where C;^ creates an electron at site i = (ix, iy) with spin 
projection a, Si=^^^ cl^cTapCip is the spin of the mo- 
bile electron, the Pauli matrices are denoted by cr, Sj is 
the localized spin at site i, (ij) denotes nearest-neighbor 
(NN) lattice sites, t is the NN-hopping amplitude for the 
electrons, J > is an antiferromagnetic (AF) coupling 
between the spins of the mobile and localized degrees 
of freedom at the same site, and J' > is a direct AF 
coupling between the localized spins in nearest neighbor 
sites. The density {n)—l — x of itinerant electrons is con- 
trolled by a chemical potential ^. Hereafter t = 1 will be 
used as the unit of energy. ^From previous phenomeno- 
logical analysis the coupling J is expected to be larger 
than t, while the Hcisenberg coupling J' is expected to 
The value of J will be here fixed to 2 



be smaller. [10 



and the coupling = 0.05, as in Ref. |7|]. To simplify the 
numerical calculations, avoiding the sign problem, the 
localized spins are assumed to be classical (with |Si| = l). 
This approximation was discussed in detail in Ref. |^] 
and it is not expected to alter qualitatively the behavior 
of the striped ground-state. The model will be studied 
using a standard MC method, details of which can be 
found in Ref. Periodic boundary conditions (PBC) 
are used. 

The spectral function A{k,uj) is defined as 



A{k,uj) 



1 



ImG{k, w). 



(2) 



where G{k,uj) is the one-particle retarded Green's func- 
tion for the electrons, k is the momentum and uj denotes 
the energy. Since we are performing a MC calculation 
on the classical spins only, because the fermions do not 
interact among themselves directly but only through the 
classical spins, the time-dependent Green's function can 
be straightforwardly calculated in real-time. The spec- 
tral functions obtained with this procedure only contain 
small controlled statistical errors. This should be con- 
trasted against other calculations that use the Maximum 
Entropy technique which is typically an uncontrolled pro- 
cedure. The density of states (DOS) N{uj) is obtained by 
adding the spectral functions for all different momenta in 
the first Brillouin zone (FBZ). 



III. SCHEMATIC PICTURE 

ARPES results in Lai.28Ndo.6Sro.i2Cu04 show a FS 
which appears to be formed by the superposition of 
two one-dimensional (ID) FS believed to be caused by 
metallic stripes in an insulating background. jlSj Simi- 
lar features are observed in underdoped LSCO p] and in 



BiaSraCaCuaOg+y (Bi2212). Q] The naive band struc- 
ture which would provide such a FS is shown in Fig.l-a, 
along the main directions in the FBZ, for the case of 
quarter-filled (metallic) equally-spaced vertical stripes, 
with a period a=4, on an AF (insulating) background. 
The overall density is (n) ~ 0.875 which corresponds to a 
hole density x = 0.125. A valence and a conduction band, 
separated by the AF gap, are associated with the insu- 
lating background. The dashed line indicates "shadow 
bands" . ||l5| In the middle of the gap there is a band that 
corresponds to the ID stripes. Notice that the mid-gap 
band is not symmetric under 7r/2 rotations. At (7r,0) 
[(0,7r)] the mid-gap band is closer to the valence [con- 
duction] band. The mid-gap band energy is given by 
Emikx,ky) — t'cosky, while the energy of the valence 
and conduction bands are given by — zti/e^ -|- 
where = —2t{coskx + cosky), and 2A is the AF gap. 
In Fig.l-a we have used t' — 0.5, t = 1 and A = 1 for 
simplicity. The metallic band is crossed by the Fermi en- 
ergy (Ep) defining the FS shown in Fig.l-b. Notice that 
the mid-gap band is flat, i.e., it does not disperse, along 
(0, tt) — (tt, tt) and (0, 0) — (tt, 0) indicating that the holes 
are static in the horizontal direction (they only can move 
vertically along the stripe) . The corresponding schematic 
density of states is shown in Fig.l-c. 
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FIG. 1. (a) Schematic band structure along the path 
(0, 0) - (0, tt) - (tt, tt) - (tt, 0) - (0, 0) - (vr, tt) in the FBZ, for 
an AF system with metallic quarter-filled vertical stripes with 
(n) = 0.875. Ef indicates the position of the Fermi energy 
and shadow bands are indicated by a dashed line. E'p indi- 
cates the position of the Fermi energy in the case of empty 
((n) = 0) stripes, (b) Fermi surface corresponding to the 
band structure shown in (a), (c) The corresponding density 
of states of Fig.l-a. 

This naive picture also shows that if the hole density 
of the stripes were 1 rather than 0.5 the system would 
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be insulating since, as shown in Fig.l-a, the Fermi en- 
ergy, indicated by Ep, would be inside a gap and all the 
stripe states would be empty. When typical mean-field 
calculations allowing for the possibility of vertical stripes 
are performed on the Hubbard model, the stable solution 
is qualitatively similar to the one shown in Fig.l-a, with 
Fermi energy E'p, indicating an insulating ground-state. 
[0 The addition of diagonal hopping terms distorts the 
bands allowing the stabilization of metallic stripes. |^ 

Although the mean-field results mentioned above are 
very instructive, it is not clear that they represent the 
true ground-state of the system. In addition, numerical 
studies are very difficult to perform directly in the Hub- 
bard model. However, we have shown that the easier-to- 
study spin-fermion model captures many of its proper- 
ties. We have also observed that the addition of neg- 
ative diagonal hopping parameters to the spin-fermion 
model with absolute values larger than 0.05i destabilize 
the striped ground-state on a square lattice. This is 
the reason why, in spite of the encouraging mean-field 
results of Ref. we will here not show data for finite 
diagonal hoppings. 
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FIG. 2. (a) Schematic band structure along the path 
(0,0) - (0,7r) - (7r,7r) - (7r,0) - (0,0) - (7r,7r) in the FBZ, 
for an AF system with a splitted mid-gap band associated to 
metallic quarter-filled vertical stripes with (n) = 0.875. Ef 
indicates the position of the Fermi energy and shadow bands 
are indicated by a dashed line. E'p indicates the position of 
the Fermi energy in the case of empty stripes, (b) The corre- 
sponding density of states. 

In order to compare the schematic results of Fig.l-a 
with those of the spin-fermion model it has to be con- 
sidered that antiferromagnetic models tend to have sym- 
metrical band structures below and above the antiferro- 
magnetic gap. Just as the conduction and valence band 
described above are caused by an AF splitting, we should 



expect that the mid-gap band will itself be splitted in the 
same way. In our schematic model the energy of the split- 
ted mid-gap band is given by = ±-y/ [t'coskyY + A'^, 
where 2A' is the gap between the two bands. The ex- 
pected band structure, using A' = 0.5 and t' — 0.7, is 
schematically shown in Fig.2-a. This indicates that, for 
vertical stripes, mid-gap spectral weight will appear close 
to the valence band around (tt, 0), while only shadow, and 
thus much weaker, spectral weight will appear around 
(0, tt). The observed mid-gap band should not disperse 
from (0,0) to (tt, 0), but it will disperse from (7r,0) to 
(tt, tt). It is also clear that very little mid-gap spectral 
weight will be observed along the diagonal direction for k 
larger than 7r/2 since the shadow bands are not intense. 
In this naive picture, the Fermi surface will be similar to 
the one shown in Fig.l-b if Ep, in Fig.2-a, is the Fermi 
energy. As before, no FS will be observed if the Fermi 
energy is E'p in Fig.2-a. The DOS, on the other hand, 
will look qualitatively different to the one displayed in 
Fig.l-c. Two mid-gap bands will appear, instead of one, 
as shown in Fig.2-b. 
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FIG. 3. (a) Schematic band structure along the path 
(0, 0) - (0, tt) - (tt, tt) - (tt, 0) - (0, 0) - (tt, tt) in the FBZ, for 
an AF system with overlapping mid-gap and valence (conduc- 
tion) bands associated to metallic vertical stripes with overall 
(n) — 0.875. Ep indicates the position of the Fermi energy 
and shadow bands are indicated by a dashed line, (b) Fermi 
surface corresponding to the band structure shown in (a), (c) 
The corresponding density of states. 

A problem with the naive schematic picture that we 
just discussed is that, as it can be observed in Fig.2-b, 
the density of states (DOS) does not have a pseudogap at 
the chemical potential which is a feature experimentally 
observed in the cuprates. In our schematic scenario, a 
pseudogap would be observed if the mid-gap states and 
the valence band overlap with each other. The resulting 
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band structure is shown in Fig.3-a. In this case, though, 
the Fermi surface would have an electron-hke contribu- 
tion coming from the valence band (doted line in Fig.3- 
b) , superimposed to a FS qualitatively similar to the one 
shown in Fig. 1-b due to the partial filling of the mid-gap 
band, denoted by the full lines in Fig.3-b. The result- 
ing DOS, with a pseudogap at the chemical potential, is 
shown in Fig.3-c. We would like to point out that al- 
though here we have presented the simplest case of band 
overlaping, it is also possible that the two bands would 
merge together in which case a single FS would be ob- 
served. 



dashed line in Fig.4-c. The ground-state is invariant un- 
der rotations in 7r/2, i.e., the features along the direction 
(0,0) — (0,7r) are identical to those along (0,0) — (tt, 0). 
The DOS showing the valence and conduction bands sep- 
arated by the AF gap is displayed in Fig.4-d. The chem- 
ical potential is in the middle of the gap indicating that 
the system is an insulator. The results shown here are to 
be expected for an AF insulator, and they demonstrate 
that the spin-fermion model with classical localized spins 
properly reproduces the main physics of such a state. 

B. (n) = 0.935 



IV. RESULTS 
A. Half-filling 
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FIG. 4. Spectral "J'Snction fo?'1'he 'undoped spin-fermion 
model on an 8 X 8 lattice as a function of lj — fi. (a) Along the 
path (0, 0) - (0, tt) - (tt, tt); (b) along the path (0, 0) - (tt, tt); 
(c) the band structure obtained from (a) and (b). The lines 
indicating the main bands are to guide the eye. Dashed lines 
indicate shadow bands. Darker points indicate higher spectral 
weight; (d) the density of states, N{u!) vs uj — fi. 

Let us consider now the similarities and differences be- 
tween the schematic results described in the previous 
section and those arising from the spin-fermion model. 
In Fig.4-a,b the spectral function along the main direc- 
tions in the FEZ is shown for the undoped (n) = 1 spin- 
fermion model on an 8 X 8 cluster. The dispersion of the 
main peaks in A(fc,a;) gives rise to the band structure 
presented in Fig.4-c, in which a valence and a conduc- 
tion band separated by the AF gap appear. The shadow 
bands associated with the AF order are indicated by the 



In the MC simulations of the spin-fermion model, holes 
added to the system tend to align forming horizontal or 
vertical stripes, as discussed in previous literature. @] 
Studying 8x8 and 12 x 12 systems it was observed that 
for X less than 1/L holes doped {L is the side of the square 
cluster) the system remains insulator. As an example, in 
Fig. 5 we present A{k, uj) along the main directions in the 
FEZ for {n) — 0.937, i.e., when 4 holes are introduced on 
the 8x8 system. 
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FIG. 5. Spectral function for the spin-fermion model with 
(n) = 0.937 on an 8 X 8 lattice, with holes aligned along the 
vertical axis forming a straight segment. The mid-gap states 
are shaded in black, (a) Along the path (0,0) — (0, tt) — (tt, tt) . 
The arrow indicates the mid-gap shadow states with nearly 
negligible weight; (b) along the path (0, 0) — (tt, 0) — (vr, tt); c) 
along the path (0,0) — (7r,7r). 
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FIG. 6. Band structure for the spin-fermion model with 
(n) = 0.937 on an 8 X 8 lattice obtained from the data in Fig. 5. 
(a) Along the path (0, 0) - (0, tt) - (vr, tt) - (tt, 0) - (0, 0) - (tt, tt) 
in the FBZ. The valence, conduction, and mid-gaps bands are 
indicated with lines to guide the eye. The dashed lines connect 
points with very small spectral weights which define shadow 
bands, (b) The density of states, N{uj) vs — /i. 

In Fig.5-a,b it can be seen that spectral weight is trans- 
ferred from the valence and conduction bands to the 
middle of the gap, creating two mid-gap bands, shaded 
in black in Fig.5-a,b, similarly to the second schematic 
model described in section III. In the results shown here, 
the 4 doped holes are aligned vertically and, as expected, 
the spectral weight associated with the holes comes from 
the valence band, and it is located around (tt, 0) in mo- 
mentum space (Fig.5-b). A redistribution of spectral 
weight, shaded in black in Fig.5-a, is also observed close 
to the conduction band around (0,7r). Only negligible 
shadow weight, indicated with an arrow in Fig.5-a, ap- 
pears close to the valence band around (0, tt). It has been 
checked that, as expected, when the holes are aligned 
horizontally the states close to the valence [conduction] 
band are located close to (0, tt) [(tt, 0)]. Also in agree- 
ment with the schematic picture, we observe that the 
lower mid-gap states disperse along (7r,0) — (tt, tt), but 
not along (0,0) — (tt, 0). A somewhat unexpected result 
though is that the low mid-gap spectral weight around 
(0, 0) is negligible, and the same occurs along the diago- 
nal direction. Thus, only states with momentum close to 
(tt, 0) and (0, tt) contribute appreciably to the mid-gap 
bands in the spin-fermion model, in contrast with the 
schematic model of the previous subsection in which all 



the momenta participated. 

The resulting band structure is shown in Fig.6-a and 
the DOS is presented in Fig.6-b. This band structure is 
reminiscent of the one described in Fig. 2 for vertical in- 
sulating ID stripes in an antiferromagnetic background. 
The chemical potential is in a gap between the lower mid- 
gap band, associated to the doped holes in the stripes, 
and the valence band. The state is insulator since the 
mid-gap bands are above the chemical potential. ARPES 
results for LSCO also detect an insulating ground- 
state at low doping. 

The symmetry under tt/2 rotations is clearly broken 
due to the vertical stripe of holes. However, independent 
MC simulations, using different random initial spin con- 
figurations, provide outputs with vertical or horizontal 
stripes indicating that the two states are degenerate. 

As more holes are added to the system, in the insulat- 
ing regime, spectral weight is transferred from the valence 
and conduction bands to the mid-gap states. 
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FIG. 7. Spectral function for the spin-fermion model with 
(n) — 0.875 on an 8 X 8 lattice. (a) Along the path 
(0, 0) — (0, tt) — (tt, tt). The arrow indicates the mid-gap shadow 
weight; (b) along the path (0,0) — (vr,0) — {n,n); (c) along 
the path (0,0) - (tt.tt). 

At (n) = 0.875, which corresponds to 8 holes on the 
8x8 size lattice, one stripe is stabilized in the ground- 
state. A very important issue is whether the system is 
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insulator, as in the early mean-field results |]T6[, or metal- 
lic as suggested by the experimental data. 

The electronic density of the stripes observed in the 
spin-fermion model is ~ 0.5, 0| in agreement with the 
experiments, but there is a slight reduction of the density 
of all the sites of the lattice indicating that, on average, 
4 holes are located at the stripe, while the other 4 are 
spread in the system. We will discuss this point in more 
detail below. Notice that this is different from the pic- 
ture emerging from experiments performed in Nd-doped 
LSCO in which two rather than one stripe would be 
expected to be observed in the situation that we are dis- 
cussing, since the holes, are assumed to be located exclu- 
sively in the quarter-filled stripes and not in the back- 
ground. This could be a drawback of this model since 
the incommensurability pattern for LSCO is similar to 
the Nd-doped material mentioned above. |J9| However, 
it is not clear that all the cuprates follow the same pat- 
tern. The spectral functions in the FBZ are shown in 
Fig. 7 and they correspond to the case in which the stripe 
is vertical. The same results, exchanging with ky, were 
obtained for the case of an horizontal stripe. In Fig.7- 
a it can be observed that the lower mid-gap states only 
have very small shadow weight, indicated with an arrow, 
around (0, tt), as expected from the schematic model, and 
the valence band remains below the chemical potential. 
Thus, it appears that there is no FS along this direction. 
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FIG. 8. (a) Eigenvalue distribution for a snapshot configu- 
ration on an 8 X 8 cluster with a vertical stripe; (b) same as 
(a) but for a 12 x 12 cluster; (c) density of states obtained 
from (a); (d) density of states obtained from (b). 

Along the (0, 0) — (tt, 0) — (tt, tt) direction, depicted in 
Fig.7-b, on the other hand, the valence band and the 
lower mid-gap band are very close to each other. On the 
8x8 system (Fig.7-b), there seems to exist a very small 
gap in between the valence and the mid-gap band which 
would indicate insulating behavior. However, we need to 
establish whether this gap is real or an artifact of the 
small lattice that we are studying. To explore this issue 



we studied the eigenvalue distribution of the Hamiltonian 
defined in Eq.(l) for several configurations of the spin 
variables (snapshots) whose ground-state had one stripe 
on 8 X 8 (Fig.8-a) and 12 x 12 (Fig.8-b) systems. We ob- 
served that the eigenvalue distribution changes very little 
for the different snapshots. On the 8x8 system a very 
small gap in the eigenvalues distribution was observed 
at the chemical potential indicated with a dashed line in 
Fig.8-a. As expected, there are 8 eigenstates forming the 
lower mid-gap states associated with the number of holes 
in the system. 
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FIG. 9. Band structure for the spin-fermion model with 
(n) = 0.875 on an 8 X 8 lattice obtained from the data in Fig. 5. 
(a) Along the path (0, 0)-(0, 7r)-(7r, 7r)-(0, 0) in the FBZ; (b) 
along the path (0, 0) — (tt, 0) — (tt, tt); (c) possible FS obtained 
from the data in Fig.7; (d) possible FS obtained by combining 
the data for vertical and horizontal stripes. The solid lines 
indicate an actual FS, while the dashed line indicates a gapped 
region or very weak FS. 

For the 12 x 12 system (Fig.8-b) the gap between the 
lower band and mid-gap states is smaller. The number 
of eigenstates per energy interval provides the density of 
states for the snapshot. An histogram of the eigenvalue 
distribution provides an estimation of the DOS. These 
histograms, using an energy interval equal to 0.04i, are 
shown in Fig.8-c for the 8x8 system, and in Fig.8-d for 
the 12 X 12 cluster. In both cases the chemical potential 
is very close to a pseudogap in the DOS located between 
the valence and lower mid-gap band, as experimentally 
observed in the cuprates. ||2l[| 

Thus, combining the results in Fig.7 and Fig. 8 we infer 
that, in this case, both the valence and the low mid-gap 
band merge together and will contribute to the Fermi 
surface. It is not possible to detect the position of the 
FS precisely due to the flatness of the overlapped bands 
close to (tt, 0). Our data are consistent with either a FS 
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at (tt - 5, 0) or at (tt, 5), where < 5 < 37r/4. 

Notice that the dispersive behavior of the mid-gap 
states is the one expected for vertical stripes in the 
schematic system, although the dispersion is small. 

A similar analysis looking at Fig.7-c shows that a FS 
along the diagonal of the FBZ i.e., from (0,0) to (tt, tt) 
does not exist. The valence band quasi-particle peak, 
which is well defined at (0,0), becomes incoherent just 
below the chemical potential for momentum (7r/2,7r/2), 
and it is almost touching the Fermi energy. The mid-gap 
states, as expected, have only shadow spectral weight. 

The band structure arising from Fig. 7 is presented in 
Fig.9-a, b and the possible FS due to the vertical stripe is 
shown in Fig.9-c. The shadowed region indicates its pos- 
sible location and the solid lines running from (77r/8,0) 
to (77r/8, tt/8) and (tt, tt/S) to {7n/8, tt/S) denote its two 
approximate extreme positions. 

In Fig.9-d we show the possible FS obtained as a super- 
position of the results for horizontal and vertical stripes. 
It is an incomplete FS which could close either around 
(0,0) or around (7r,7r) (solid line in Fig.9-d). We could 
not determine if the surface is close or open but away 
from (tt, 0) and (0, tt) the spectral weight is weak and in- 
coherent. The possible extreme positions of this weaker 
FS (or gap) is indicated with dashed lines in Fig.9-d. 
The FS closing around (tt, tt) has been observed in un- 
derdoped LSCO @ and Bi2212 The observation 

of a FS closing around (0, 0) is controversial. According 
to Ref. it ap£ears in Bi2212 at 32.3 eV photon en- 
ergy, while Ref. |22[| interprets similar data as indicating 
no changes in the shape of the FS. 

Our results are also in agreement with the LSCO data 
showing a weak or no FS along the diagonal in the un- 
derdoped regime. |^ ARPES data seem to indicate a FS 
along the diagonal direction in Bi2212 but we be- 
lieve that inverse ARPES data would be necessary to rule 
out the existence of a gap. In our results the existence 
of the gap becomes clear after observing the behavior of 
the spectral weight in the valence band and the mid-gap 
band above ^. 



D. (n) = 0.80 

Next we will study how the band structure evolves as 
the system is doped further away from half-filling. At 
^ 20% doping it is possible to determine that the FS 
closes around (0,0). The symmetry under rotations is 
still spontaneously broken, but the potential barrier sep- 
arating vertical and horizontal stripes appears to be lower 
than for lighter doping. As a result, some features due to 
horizontal configurations appear in the spectral functions 
presented here which corresponds to a ground-state with 
mostly vertical charge inhomogeneities. 

The spectral functions along the main directions in the 
FBZ are shown if Fig.lO-a, b and c, and the correspond- 
ing band structure is in Fig.lO-d and e. As in the previous 



case, a gap rather than a FS is observed close to (0, tt) 
(Fig.lO-a and d). |Q Close to (tt, 0) the valence and mid- 
gap bands overlap with each other where the FS seems to 
be located (Fig.lO-b and e). The FS along the diagonal 
is still inexistent or very weak since the spectral weight 
at (7r/2,7r/2) is incoherent (Fig.lO-c). 
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FIG. 10. Spectral functions for the spin-fermion model 
with (n) = 0.81 on an 8 x 8 lattice. (a) Along the 
path (0,0) - (0,7r) - (tt.tt) in the FBZ; (b) along the path 
(0,0) - (7r,0) - (7r,7r); (c) along the path (0,0) - (7r,7r); (d) 
band structure along the path (0, 0) — (0, tt) — (vr, n) — (0, 0); 
(e) along the path (0,0) — (7r,0) — (7r,7r); (f) FS obtained by 
combining the data for vertical and horizontal stripes. The 
dashed line indicates a very weak FS; (g) The density of states. 

The shape of the FS for the combined ground-state 
(with vertical and horizontal incommensurability) is 
shown in Fig.lO-f. This FS was obtained by analyzing 
A{k,uj) for all the available momenta in the FBZ. In 
Fig.lO-g we present the DOS and it can be seen that 
the chemical potential is in a pseudogap. 



E. (n) = 0.75 

When 16 holes are introduced in the system two par- 
allel stripes are stabilized. The system is metallic and 
the FS closes around (0,0). The spectral functions along 
(0, 0) — (0, tt) — (tt, tt) are shown in Fig.ll-a. It is clear that 
the valence band does not cross the chemical potential, 
and only shadow weight appears in the lower mid-gap 
band. Along (0,0) — (7r,0) — (tt, tt), presented in Fig. 11- 
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b, the valence and low gap bands overlap and cross the 
chemical potential at A: ~ (37r/4,0). The FS in the diag- 
onal direction is weak because the spectral weight that 
crosses to — fj, is incoherent as it can be seen in Fig.ll-c. 
The band structure corresponding to the case in which 
the stripes are vertical is shown in Fig.ll-d and e. The 
superposition of the two ground-states (vertical and hori- 
zontal stripes) will give raise to a FS which closes around 
(0,0) as shown in Fig.ll-f. Clearly the chemical potential 
is located in a pseudogap in the density of states, as it 
can be seen in Fig.ll-g. 
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FIG. 11. Spectral functions for the spin-fermion model 
with (n) = 0.75 on an 8 x 8 lattice. (a) Along the 
path (0,0) - (0,7r) - (tt.tt) in the FBZ; (b) along the path 
(0,0) - (7r,0) - (7r,7r); (c) along the path (0,0) - (7r,7r); (d) 
band structure along the path (0, 0) — (0, tt) — (vr,7r) — (0, 0); 
(e)along the path (0,0) - (vr, 0) - (7r,7r) (f) FS obtained by 
combining the data for vertical and horizontal stripes. The 
dashed line indicates a very weak FS; (g) The density of states. 



F. (n) = 0.625 

With increasing doping the area of the FS is reduced. 
For (n) = 0.625 the symmetry under rotations is still bro- 
ken. The magnetic and charge incommensurability still 
form vertical or horizontal patterns. The spectral func- 
tions for the vertical pattern are shown in Fig.l2-a, b 
and c. In Fig.l2-a it can be seen that no FS is observed 
along (0,0) - (0,7r) - (7r,7r). Fig.l2-b shows that the 
mid-gap and valence band which had merged together 
cross /i below (37r/4,0). In the diagonal direction just 



below (7r/2, 7r/2) incoherent weight crosses the chemical 
potential determining a weak FS (Fig.l2-c). The corre- 
sponding band structures are presented in Fig.l2-d and e 
and the FS is shown in Fig.l2-f. The chemical potential 
is in a pseudogap as shown in Fig. 12-g. 
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FIG. 12. Spectral functions for the spin-fermion model 
with (n) = 0.625 on an 8 x 8 lattice. (a) Along the 
path (0,0) - (0,7r) - (tt.tt) in the FBZ; (b) along the path 
(0,0) - (7r,0) - (7r,7r); (c) along the path (0,0) - (7r,7r); (d) 
band structure along the path (0, 0) — (0, vr) — (vr, vr) — (0, 0); 
(e)along the path (0,0) - (vr.O) - (tt.tt) (f) FS obtained by 
combining the data for vertical and horizontal stripes. The 
dashed line indicates a very weak FS; (g) The density of states. 



G. Analysis 

The previous results indicate that as the electronic 
density varies from 1.0 to 0.625 a clear change in the 
band structure of the system occurs. Spectral weight for 
the momenta close to the chemical potential is transfered 
from the valence band to the lower mid-gap band. As a 
result, the chemical potential always appears in a pseu- 
dogap between the valence and the lower mid-gap bands. 
The bottom of the valence band, located at fc = (0,0), 
remains at an approximately constant distance from the 
chemical potential. This is very different from the mean- 
field studies 1^ in which a rigid band structure is ob- 
tained and the chemical potential approaches the bot- 
tom of the valence band as holes are added. We believe 
that the main difference between our numerical results 
and mean-field is due to the fact that in the mean-field 
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analysis the proposed striped background is always an 
almost perfect antiferromagnet, which physically is not 
accurate. Another important characteristic of the evolu- 
tion of the band structure with doping is that while the 
gap between the valence and lower mid-gap band closes 
becoming a pseudogap, the gap between the valence and 
conduction band still exists even at 40% doping. 

Another characteristic of the band structure in the 
metallic case is that for vertical incommensurability, the 
quasiparticle at (0, tt) has lower energy than at (tt, 0) as 
it can be seen in Fig.9-a,b, and in parts (d) and (e) 
of Figs. 10, 11 and 12. If we consider a 2D system of 
free electrons with horizontal hopping tx and vertical 
hopping tj,, the energy as a function of momentum is 
given by ek = —2{txCos{kx) + tyCos{ky)). It is clear that 
e(7r, 0) > e(0,7r), if tx > ty. This may indicate that, in 
the spin-fermion model, it is easier for the electrons to 
move in the direction perpendicular to the stripes rather 
than along them. 

V. THE PSEUDOGAP 

The results presented in the previous section show 
that, upon hole doping, a pseudogap between the valence 
and the lower mid-gap stripe-induced band develops at 
the chemical potential in the density of states (see Fig.8- 
c,d, and part (g) of Figs. 10, 11 and 12). This is a fea- 
ture that has been observed experimentally both in the 
cuprates and in the manganites , and it deserves 
to be understood in terms of the spin-fcrmion model. 

The formation of a pseudogap at the chemical potential 
in the manganites has been explained in previous stud- 
ies as arising from the coexistence of hole-poor antiferro- 
magnetic insulating regions and hole-rich ferromagnetic 
clusters. ^ In Rcf. |g5| it was concluded that pseudo- 
gap behavior should be observed in any compound which 
is in a mixed-phase regime. Although the spin-fermion 
model does not present phase separation for the value 
of J studied in the present paper, it is clear that the 
"striped" ground state has mixed-phase characteristics, 
since the stripes are richer in holes than the background. 
A schematic picture of the hole distribution in the di- 
rection perpendicular to the stripes is shown in Fig.l3-a. 
The density of holes x in the direction perpendicular to 
the stripes is clearly non-uniform, and it is schematically 
shown in Fig.l3-b. We can represent this result using 
an "effective potential" for the holes which will be more 
negative at the stripes, as shown in Fig.l3-c. The states 
inside the wells generate the mid-gap bands along the di- 
rection parallel to the stripes in the BZ, (7r,0) — (tt, tt) 
for vertical stripes. Clear quasi-particle peaks and maxi- 
mum dispersion are expected along this direction since it 
is easier for the holes to move along the stripes. Along the 
direction perpendicular to the stripes, from (0, 0) to (tt, 0) 
for vertical stripes, the holes move through tunneling be- 
tween the potential wells or by fluctuations away from 



the stripe, effectively providing a width d to such stripe. 
This motion will generate a mid-gap band, schematically 
shown in Fig.l3-c. The electrons, on the other hand, will 
move mostly across the stripes, as discussed in the pre- 
vious section and thus, a dispersive valence band would 
be expected. As a result of the formation of the mid-gap 
bands a pseudogap develops at the chemical potential in 
the DOS, as schematically shown in Fig. 13-d. 




Band 

FIG. 13. Schematic explanation of the pseudogap forma- 
tion upon hole doping in the spin-fermion model. In (a) a 
typical distribution of the AF background and stripes (S) 
reach in holes along the direction perpendicular to the stripes 
is shown. In (b), the corresponding hole density x is sketched, 
showing that the holes are mostly located in the stripes. In 
(c), the effective potential felt by the holes is presented. The 
thick line indicates the mid-gap band populated with holes 
which develops. The resulting DOS with the chemical poten- 
tial in the pseudogap is sketched in (d). 

An important difference between the pseudogap in the 
spin-fermion model and the one observed in models for 
manganites is that in the latter the pseudogap is isotropic 
in momentum space. In the spin-fermion model, on the 
other hand, the pseudogap is clearly observed close to 
(tt, 0), while a clear gap appears close to (0, tt) (for verti- 
cal stripes) since the lower mid-gap band has only shadow 
weight here and the valence band has a lower energy than 
at (tTjO). This difference should be attributed to the 
one-dimcnsional nature of the stripes in the spin-fermion 
model. The incoherence of the spectral weight along the 
diagonal is, as mentioned above, due to the fact that di- 
agonal hopping of electrons or holes is not particularly 
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favored. 



VII. ACKNOWLEDGEMENTS 



VI. CONCLUSIONS 

Summarizing, we have studied the spectral functions 
of the spin-fermion model, which has a ground-state in 
which added holes tend to form vertical and horizontal 
stripes, using unbiased numerical techniques. We have 
observed that doped holes contribute to the formation of 
mid-gap bands by modifying the valence and conduction 
bands associated to the insulator. The mid-gap bands 
arise from a dynamically generated effective potential 
which produces a non-uniform distribution of charge in 
the ground-state. In the metallic regime the lower mid- 
gap and valence band overlap with each other, giving rise 
to a pseudogap in the density of states at the chemical po- 
tential, as observed in experiments for the cuprates. A 
pseudogap arising from inhomogeneities in the ground- 
state has also been observed in previous investigations of 
models for manganites. 

The ground-state of the spin-fermion model appears to 
change from an insulator to a conductor with an incom- 
plete FS when a single stripe gets stabilized, i.e., for a 
hole density x — \/L. In the conductor, electrons hop 
more easily in the direction perpendicular to the stripes. 
Thus, though the striped state in the spin-fermion model 
is not insulating, it appears that most of the electronic 
hopping occurs in the direction perpendicular to the 
stripes rather than along the stripes. Close to 20% dop- 
ing the FS clearly closes around (0,0). These features are 
qualitatively similar to those observed in ARPES exper- 
iments. The band structure is not rigid, in disagreement 
with mean-field results, and it smoothly changes from 
insulating to conductor. Another interesting feature is 
the observation of well defined quasi-particle peaks in 
A{k,uj) which become incoherent as the chemical poten- 
tial is approached indicating that the FS is very weak or 
it does not exist along the diagonal direction. Also close 
to (0, tt) [(tt, 0)] in the case of vertical [horizontal] stripes 
the chemical potential is in a gap. These features are 
related to the one-dimensionality of the stripes. A sim- 
ilar behavior of the valence and mid-gap band has been 
observed using a cluster perturbation theory for the t — J 
model. Ill 

Although part of our results disagree with some of the 
theoretical models for stripes in the cuprates, we believe 
that this numerical study provides unbiased information 
on the dynamical properties of charge inhomogeneous 
ground states which will be relevant to understand the 
behavior of materials whose inhomogeneous properties 
are just beginning to be unveiled experimentally. 
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